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ABSTRACT 

We examine the case for Quark-Novae (QNe) as possible sources for the reionization and early 
metal enrichment of the universe. Quark-Novae are predicted to arise from the explosive collapse 
(and conversion) of sufRciently massive neutron stars into quark stars. A Quark-Nova (QN) can occur 
over a range of time scales following the supernova event. For QNe that arise days to weeks after 
the supernovae, we show that dual-shock that arises as the QN ejecta encounter the supernova ejecta 
can produce enough photons to reionize hydrogen in most of the Inter-Galactic medium (IGM) by 
z ~ 6. Such events can explain the large optical depth Tg ~ 0.1 as measured by WMAP, if the 
clumping factor, C, of the material being ionized is smaller than 10. We suggest a way in which a 
normal initial mass function (IMF) for the oldest stars can be reconciled with a large optical depth 
as well as the mean metallicity of the early IGM post reionization. We find that QN also make a 
contribution to r-process element abundances for atomic numbers A > 130. We predict that the main 
cosmological signatures of Quark-Novae are the gamma-ray bursts that announce their birth. These 
will be clustered at redshifts in the range z ~ 7-8 in our model. 
Subject headings: Reionization - Quark-Nova - r-process - Gamma- Ray bursts 



1. INTRODUCTION 

The end of the cosmic "dark age" began with the pro- 
duction of reionizing UV radiation from the first lumi- 
nous objects. Several independent lines of evidence, in- 
cluding recent WMAP results suggest that beginning at 
redshifts of about z 17, an inhomogenous and possi- 
bly non-monotonic decrease in the fraction of neutral H 
and He in the IGM took place, and that this process was 
largely completed by z '~ 6. The nature of the first ion- 
izing sources is still mysterious - possible candidates for 
"early" sources (i.e., well before z ~ 6), include metal- 
free Population III (Pop III) stars, high-z dwarf galaxies, 
intermediate-mass black holes or sterile neutrinos. Cur- 
rent models favour luminous objects that formed from 
collapsed gas haloes of mass M ~ IO^Mq. 

Although quasars and active galactic nuclei (AGNs) 
are efficient emitters of UV photons, results from the 
Sloan Digital Sky Survey (SDSS) imply that they con- 
tribute much less to the ionizing radiation than star- 
forming galaxies at high redshift (Fan et al. 2006). The 
first stars, owing to their metal-free composition (Z=0) 
which makes them a copious source of ionizing photons, 
are increasingly viewed as prime candidates for early 
reionization sources. It has been argued that a top-heavy 
IMF (40-lOOAfo) may be required to explain the large 
WMAP optical depth and the ongoing overlap of HII 
regions at lower redshifts z < 9 suggested by IGM tem- 
peratures inferred from the Lya forest (Theuns et al. 
2002). Pop III stars derived from a top-heavy IMF can 
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emit as many as 10^ ionizing photons per baryon dur- 
ing the star's life (Tumlinson&ShuU 2000; Bromm et al. 
2001; Schaerer 2002). A population of very massive stars 
(VMS), that have been linked to pair-instability SNe, 
was advocated by Qian & Wasserburg (2002) to explain 
r-process element abundances as well as that of C, O and 
Si in some extremely metal-poor stars, while other works 
(Tumlinson et al. 2004; Daigne et al. 2004; Venkatesan 
& Truran 2003; Fang&Cen 2004) argue that an ordinary 
Type II-SN IMF can explain these abundances just as 
well while also providing the large optical depth Te mea- 
sured by WMAP. Furthermore, the escape efficiency of 
UV photons required for early reionization is smaller for 
a top-heavy IMF than for the VMS, making the former 
better candidates for reionization sources. 

While the high production of ionizing photons from 
a top-heavy Z—0 stellar population is desirable, their 
short lifetimes (few million years in the absence of sig- 
nificant mass loss) imply that successive generations do 
not remain metal-free. This creates some major diffi- 
culties. The issue here is that the IGM is rapidly pol- 
luted in metals ^. This "negative feedback" mechanism 
(e.g. Sokasian et al. 2004; Ricotti & Ostriker 2004) is a 
problem - it might result in the reneutralization of the 
IGM, calling for another sub-epoch of reionization from 
a different stellar population in a metal-enriched IGM. 
This may necessitate double or even multiple reioniza- 
tion (Wyithe & Loeb 2003). It may be possible to extend 
the Pop III lifetime by "hiding" most heavy elements in 
the collapse to black holes, with added benefits of ad- 

* Analysis of quasar metal absorption lines have shown that the 
IGM is enriched in heavy elements from a level of ~ 10~^ Zq at 
z ~ 5.3 to ~ 10"'^ Zq at 2 < 4.3 where Zq is the solar metallicity 
(Songaila 2001). That such an abundance of heavy elements is 
observed commonly and uniformly in the low-density IGM where 
no bright galaxies are found may be due to outflows from small 
protogalaxies (Madau et al. 2001; Mori et al. 2002), but prompt 
enrichment by Pop III stars provides a more natural solution for 
the distribution of metagalactic heavy elements. 
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ditional accretion generated X-ray photons, but this ap- 
pears to require extreme black hole formation rates and 
smaller than normal metal yields from Pop II stars (see 
Ricotti&Ostriker 2004 for a discussion). 

In this paper, we argue that explosions that result 
from the conversion of neutron stars (born from progen- 
itors with mass in the 2^Mq < M < 40A'/q range) to 
quark stars - known as Quark-Novae (QNe; Ouyed et al. 
2002) - may be important sources of reionizing UV pho- 
tons. This mechanism has the advantage of not leading 
to heavy chemical pollution of the IGM or ISM, thereby 
alleviating the requirement of producing several gener- 
ations of chemically pristine massive stars early on in 
cosmic star formation history. The first QNe in the uni- 
verse can occur very soon after Pop III stars end their 
lives as Supernovae (SN) - in fact the time delay between 
the SN and the subsequent QN is one of the main param- 
eters of QN theory (Ouyed et al. 2007). The ioinizing 
photons are produced when QN shock waves overtake 
their progenitor SN shocks (Leahy&Ouyed 2008). Co- 
pious energy releases can occur in this way if the QN 
delay is days to months after the initial SN. QNe can 
provide enough UV photons to explain the large optical 
depth measured by WMAP. The contribution to reioniza- 
tion from the main-sequence lifetime of the first stars is 
then reduced, which eliminates the need for a top-heavy 
IMF. The "negative feedback" problem of the increasing 
metallicity is mitigated by effectively decoupling metal- 
enrichment from reionization in the QN picture, without 
recourse to a large formation rate of black holes. QNe 
also provide a local and prompt r-process enrichment for 
elements above A ^ 130, which can be related to the ele- 
ment abundance patterns in extremely metal-poor stars 
once chemical evolution studies are performed. 

The properties of QN summarized above are the basis 
of our proposed modification of reionization and chem- 
ical enrichment wherein: (i) the most massive stars 
(M > 4OM0) in a conventional or slightly top-heavy 
IMF collapse to black-holes with a possible (small) con- 
tribution to reionization from accretion but no contri- 
bution to metallicity while (ii) the reionization is driven 
by intermediate-mass Pop III stars, whose higher mass 
members (25Mq < M < 40Mq) end up as QNe, pro- 
viding the bulk of reionzing photons and enriching their 
environment in elements beyond A ~ 130, and (iii) low 
mass Pop III stars (SMq < M < 25Mq) end up as type 
II SNe. A long-lived and metal-poor population of low- 
mass stars begins to emerge at the end of the reionization 
epoch (Greif et al. (2008)). In our scenario, the dying 
out of the first heavy stars coincides with a peak in the 
QN rate and therefore a peak in ionizing radiation. 

In this paper, we first outline the basic observational 
constraints on reionization models (§2), followed by a 
description of QNe features in §3 (the explosion, the 
compact remnant and the ejecta), and then derive the 
properties of reionization by QNe (§4) as well as early r- 
process enrichment (§5). Our discussion and conclusions 
(§6) suggest observational tests of this scenario. 

2. OBSERVATIONAL CONSTRAINTS ON REIONIZATION 

We first examine the constraints provided by the 
Gunn-Peterson test and the latest WMAP results. Since 
even a tiny fraction of neutral H (cchi ~ 10"'^) in the 
IGM is sufficient to extinguish the transmitted flux blue- 



wards of the Lya emission line from distant Quasars, the 
onset of the Gunn-Peterson trough (Gunn & Peterson 
1965) in recent observations of Quasar emission spectra 
(Schneider, Schmidt, & Gunn 1991; see Djorgovski et al 
2001 and references therein) and the increased variance 
in the Gunn-Peterson optical depth at 2; ~ 6 shows that 
the IGM was transiting rapidly to a fully reionized state 
by then. This observation, although not a deep probe 
of reionization history, places a floor on the emissivity 
of UV photons required to keep up with recombination 
during the reionization epoch (Miralda-Escude, Haehnet 
& Rees 2000): 

*,„..„ = 10".'.-.Mpe-3(g)(i±i)=(0^)\ 

(1) 

where h = i7o/(100 km s ^ Mpc ^) is the scaled 
Hubble constant, and C is the clumping factor of the 
IGM related to smaU-scale (~ kpc) gas inhomogeneities 
caused by structure formation (Shapiro&Giroux 1987). 
We adopt a ACDM cosmology with h = 0.72, fioM = 0.26 
and r^A = 0.74, ilb/i^ ~ 0.023 which are consistent with 
the 5-year WMAP data (Dunkley et al. 2009). 

The WMAP results provide important constraints on 
the Thompson optical depth at the time of reioinization. 
Cosmic Microwave Background (CMB) photons in the 
reionized epoch are Thomson rescattercd by free elec- 
trons. This damps the primary temperature anisotropy 
of the CMB on scales smaller than the horizon size 
at rescattering, and generates additonal polarization on 
large scales (tens of degrees). The WMAP mission's mea- 
surements of the CMB polarization spectrum at low I 
have quantified this effect through an integrated Thom- 
son optical depth Tg. Recent analysis of 5-year WMAP 
data constrains it to be Tg = 0.087 ± 0.017 for the best- 
fit ACDM cosmology (Dunkley et al. 2009); this is a 
much tightened value compared to the 3-year WMAP 
value Te = 0.089 ± 0.030 (Spergel et al. (2007)). These 
results rule out prompt reionization at redshift 2 ~ 6 at 
3.5(7 and taken together with the varying Gunn-Peterson 
optical depth at z ^ 6, argue for an extended period of 
reionization stretching back to z ~ 20. The smaller mass 
fraction in collapsed halos at such high z (e.g. Iliev et 
al. 2007 and references therein) implies that a high emis- 
sivity of ionizing photon per unit mass is required early 
in the reionization epoch. This may be achieved with 
high formation efficiency of metal-free massive Pop III 
stars, but as we show in sections 3 and 4, our estimates 
for QN formation rates, the QN yield in heavy elements 
and UV photons imply that QNe allow for a more con- 
ventional star formation rate for Pop III stars as a viable 
alternative. 

3. THE QUARK-NOVA 

If a neutron star undergoes a deconfinement phase 

transition to quark matter in its interior, it can im- 
dergo core collapse, resulting in mass ejection of the 
outer crust (Takahara & Sato 1986; Gentile et al. 1993; 
Fryer & Woosley 1998;) and the formation of a stable 
strange quark star (a star made of up, down and strange 
quarks; e.g. Itoh 1970; Bodmer 1971; Witten 1984; Al- 
cock et al. 1986, Bombaci et al. 2004). Ouyed et al. 
(2002) suggested that for sufficiently massive progeni- 
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tors, the dcconfincd cores physically separate from the 
overlaying hadronic layers giving rise to powerful no- 
vae. They termed this the Quark-Nova scenario (see also 
Keranen&Ouyed 2003). 

As outlined in Keranen et al. 2005, the initial state for 
the QN is that of a deleptonized neutron star with a (u,d) 
core. In the quark-nova (QN) picture, the (u,d) quark 
core of a hybrid star that undergoes the phase transition 
to the (u,d,s) quark phase, shrinks in a spherically sym- 
metric fashion to a stable, more compact strange mat- 
ter configuration faster than the overlaying material (the 
neutron-rich hadronic envelope) can respond, leading to 
an effective core collapse. The two-step process, neutron 
to (u,d), then (u,d) to (u,d,s) is crucial. In this scenario, 
the neutrinos come from weak reactions at the edge of 
the (u,d) core and can leak out easily into the surround- 
ing cooler and deleptonized envelope where they can de- 
posit energy. This is significantly different from phase 
conversion in a proto-neutron star stage where neutrino 
transport is slower (of the order of seconds) because of 
the hot and lepton-rich matter. 

A complete dynamical treatment of the QN, including 
neutrino transport and full stellar evolution is only just 
beginning to be explored (Sagert et al. 2009) and pre- 
liminary results support a strong secondary explosion. 
However, we can make reliable estimates by adopting 
the viewpoint of numerous previous authors (Horvath & 
Benvenuto 1988; Lugones et al. 1994; Anand et al. 1997; 
Ouyed et al. 2002; Keranen et al. 2005; Jaikumar et al. 
2007; Sagert et al. 2009) who assume a rapid adiabatic 
collapse on gravitational free-fall timescales (~ 10~*s) 
along with a conversion of the core to (u, d, s) quark 
matter on weak interaction timescales (^ 10~*s). The 
gravitational potential energy released in the collapse 
is converted partly into latent heat of the (presumably) 
first-order phase transition and partly into outward prop- 
agating shock waves from the supersonic motion of the 
(u, d, s) conversion front ^ . The temperature of the quark 
core thus formed rises quickly to (10-20) MeV since the 
collapse is adiabatic rather than isothermal (Gentile et 
al. 1993). 

3.1. The quark-nova compact remnant 

We assume that hot quark matter in the color-flavor- 
locked (CFL) phase is the true ground state of matter 
at high density. This is a superconducting phase that is 
energetically favored at extremely high densities and low 
temperatures (CFL; Alford et al. 1999). In this phase 
u, d, and s quarks pair, forming a quark condensate (a 
superfluid) that is antisymmetric in color and flavor in- 
dices. This state is reached by the QN compact remnant 
as it cools below a few tens of MeV (see Appendix A for 
a disciission on CFL and contamination of the universe) . 

Mannarelli et al. (2008) argue that CFL stars are un- 
likely to constitute a significant number of puslars since 
the r-mode is undamped for frequencies above IHz, due 
to the weak damping effect of mutual friction. How- 
ever, this conclusion is derived for low temperatures 
T < O.OlMeV, while the temperature for the quark-nova 

^ More recent work, assuming realistic quark matter equations 
of state, argues for strong deflagration (Drago et al. 2007), but 
this too may be preceded by a compression shock (Lugones et al. 
1994) in the hadronic phase. 



is tens of MeV. As shown in Jaikumar et al. (2008), bulk 
viscosity from Kaons in the CFL phase is effective in 
damping out the r-mode at T > IMeV, allowing the star 
to spin rapidly. Thus, the r-mode instability argument 
does not contradict the quark-nova mechanism (and the 
related CFL firebaU; e.g. Vogt et al. 2004). However, 
it is clear that following the cooling phase of the quark- 
nova, the temperature will fall below 0.01 MeV. While 
mutual friction is too weak in this regime, other sources 
of viscosity may be important, such as shear viscosity 
from electromagnetic processes (it was shown in Manuel 
ct al. (2005) that shear viscosity from phonons, though 
large, is inconsequential since it violates the basic hydro- 
dynamic assumptions at such low temperatures). There- 
fore, we cannot comment on the ultimate fate of the CFL 
star below 0.01 MeV, but the quark-nova mechanism, 
which operates at much larger temperature, is not in ob- 
vious conflict with the r-mode instability. 

3.2. The quark-nova ejecta 

The gravitational potential energy released (plus latent 
heat of phase transition) during this event is converted 
partly into internal energy and partly into outward prop- 
agating shock waves which impart kinetic energy to the 
material that eventually forms the ejecta. Unlike Super- 
novae, neutrino- driven mass ejection in Quark-novae is 
not feasible, as neutrinos are trapped inside a hot and 
dense expanding quark core, once it grows to more than 
^2 km (Keranen et al. 2005). Mass ejection due to 
core bounce is also unlikely unless the quark core is very 
compact (l-2)km. A more promising alternative is mass 
ejection from an expanding thermal fireball that is a di- 
rect consequence of dense, hot quark matter in the color- 
flavor-locked (CFL) phase (Jaikumar et al. 2002; Vogt 
et al. 2004; Ouyed et al. 2005). Depending on the con- 
version efficiency of photon energy to kinetic energy of 
the ejecta, up to 1O~^M0 of neutron-rich material can 
be ejected at nearly relativistic speeds (Ouyed&Leahy 
2009). Thus, QNe ejecta decompress material from the 
outer layers of exploding neutron stars. This is impor- 
tant for the operation of the r-process (Jaikumar et al. 

2007) . 

3.3. Dual-shock quark-novae 

If a QN occurs a few days to weeks (Yasutake et al. 
2005; Stafl[ et al. 2006) after the Supernova, the QN 
ejecta propagating into the Supernova envelope rapidly 
sweep up enough mass to become sub-relativistic. The 
collision of the two ejecta sets up a blast wave that prop- 
agates outward and re-energizes the Supernova ejecta. 
Part of the reshocked ejecta's energy is then radiated 
away in photons. This model for a "dual-shock" QN 
(hereafter dsQN) has been successfully applied to the 
observed light-curve of the most energetic Supernovae 
(eq. SN2006gy, SN2005ap, SN2005gj; see Leahy&Ouyed 

2008) . It is this radiation coming from the reshocked 
ejecta which we wish to evaluate as a potential reion- 
izing source, keeping in mind that such dsQNe may be 
fairly frequent (discussed below). 

4. QUARK-NOVAE AND REIONIZATION 

4.1. The Ionization Ejfi,ciency 

We estimate the frequency of dsQNe, denoted by 
/dsQN, by assuming a conventional IMF for massive stars 
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(Scalo 1986), an average Galaxy density ng{z)—ng{Q)(l + 
z)'^, and an average Type II supernova rate fsN{z) 
per year per Galaxy. For purpose of estimation, we 
take (1 + 2;) ~ 10 as typical of the reionization epoch, 
ng(0)=0.02/(Mpc)3 and /sAr(-z) - 1 yr"! Galaxy-^. We 
find /dsQN = 0.1xQN/siv"s(0)(10)3 ~ 0.2 yr"! Mpc'^ 
where xqn 0.1 is the fraction of stars in the 25-4OM0 
range relative to all supernovae (8-4OM0) in the Scalo 
IMF. These stars are expected to undergo a Quark-Nova 
^ with delay time between the supernova and the quark- 
nova decreasing with increasing mass (we assume conver- 
sion due to increased fall-back for more massive stars; see 
also Appendix B). The additional pre- factor of 0.1 is the 
fraction of quark-novae that can support an efhcient dual 
shock. As argued in Leahy & Ouyed (2008), to be consis- 
tent with the fraction of superluminous supernovae such 
as SN2006gy, the fraction of quark-novae which happen 
fast enough for the dual shock to be effective is typically 
about 0.1. Galaxy mergers or a longer ionization history 
can easily increase /dsQN ~ 10 yr~^ Mpc"'^. We choose 
/dsQN ~ 1 yr~^ Mpc~^ as an average frequency of QN 
events during the reionization epoch. 

To determine the number of ionizing photons, we note 
that the time-dependent luminosity of a dsQN is given 
by (Leahy&Ouyed 2008) 



isN(0 = CvATc 



core "-cjcc 



47ri?phot.(i)' 



,dD{t) 
dt 



(2) 



where Cv ~ (3/2) fee is the specific heat of the ejecta, 
ATcore ~ Tcore is the core temperature of the ejecta, 
^ejec IS the number density of the ejecta, Rphot.{t) is 
the photospheric radius and D{t) is the photon diffusion 
length. This yields an integrated luminosity of lO^^ergs, 
in agreement with observations of SN2006gy, for the pa- 
rameter choice i?phot.(i) ~ 3 x lO^^km ^. Furthermore, 
during time when the ejecta are optically thick (tens of 
days for the most energetic Supernovae), the luminosity 
can be approximated by LsN(i)=cr47ri?phot(0^^^ with 
T Ki lO'^K. This is close to the observed spectral peak of 
SN2006gy, in the first tens of days( Smith et al. (2007)). 
For metal-free Pop III progenitors however, we expect 
the peak to he in the UV region (Tpk « (2 - 4) x IQ-^K). 
The number of ionizing photons is then 



No 
2.4 



dx 

13.6/Tpk(oV) & 



1 



(3) 



where A^o=0.244(47r -Rphot/3)?'pk is the total number of 
radiated photons. We find A^o=2 x 10^^-2 x 10^^ .^^i^jig 
A^ion ~ 3 X 10^^-4 X 10^^ for the lower and upper limit 
of Tpk respectively. Since we have neglected the ion- 
izing photon flux after the supernova ejecta becomes 
transparent, our estimate is a very conservative one. 

6 Staff et al. (2006) estimated that 1 out of a 1000 neutron 
stars could undergo the deconfinement transition during a Hubble 
time from spin-down alone. Mass fallback is more efficient at rapid 
conversion and upto 1 in 10 neutron stars with progenitors above 
25M0 can become quark stars. Although this rate may seem high, 
Leahy & Ouyed (2007) have argued that this is consistent with the 
inferred birth rate of AXPs and SGRs (Gil&Heyl 2007), which are 
quark stars in their model. 

i?phot. (*)=^0+'i^SJvi with i?o being the radius of the progenitor 
star and vs n being the speed of the shocked material. 



The exact number of ionizing photons is very sensitive 
to Tpk and i?phot(0; so we choose i?phot(i)=3x lO^^cm, 
Tpk=3 X lO^K, giving A^ion=5 x 10^° as a reasonable esti- 
mate (see also Appendix C). Then, the corresponding 
emissivity of UV ionizing photons is conveniently ex- 
pressed as 



A^i, 



ioni..QNc ' 



1.5 X lO'^^ s"^ Mpc-3 

/dsQN 

1 yr-1 Mpc" 



(4) 



N, 



5 X 10^° photons 



which exceeds the floor on the emissivity of ionizing 
photons from Equation ([1]), even for the lower limit on 



A^ln 



3 X 10 



59 



It follows that the total number of ionizing photons 
•per baryon produced in dsQNe during the epoch St^c (be- 
tween 2; ~ 17 and z ~ 6) is 



5 X 106") 



/dsQN 



1 Gyr y VI Gpc 



lyr-i Mpc 
3 



-3 



(5) 



1079 



where i?univ. IGpc is the comoving radius of the uni- 
verse at the reionization epoch ((1 -|-z) ~ 10) and Auniv. 
is the total number of baryons of the universe, estimated 
to be ~ 10^^ baryons. Taking into account the uncertain- 
ties in Supernova rates in the distant past, we flnd from 
Equation (O that QNe generate about 0.1-10 photons 
per baryon in the universe. A slightly higher (lower) 
Tpk or -Rphot can substantially increase (decrease) this 
percentage. However, Equation (O suggests the deflnite 
possibility that QNe can be an important source of rcion- 
izing photons. 

4.2. The Thomson Optical Depth 

The Thomson optical depth is the observational link 
with the ionization history of the universe. It is defined 
by the relation 



ne{z)(JTcdt = / f{z)n-ii{z)aTC 



dz 



{l + z)H{z) ' 

where ne{z) is the co-moving electron density and <7t = 
6.6524 X 10^^^ cm^ is the Thomson cross-section. In the 
second part of Equation ([6]), f{z) is the hydrogen ioniza- 
tion fraction, nn is the number density of Hydrogen, and 

1 /2 

H{z)=Ho (noM(l + z f + Qa) ' is the Hubble parame- 
ter in terms of the cosmological redshift in a fiat universe 
with a cosmological constant. We calculate f{z) for QN, 
and then show how observational constraints on the op- 
tical depth provide constraints on QN associated with 
Pop III stars. 

Each dsQN will set up an ionization front that prop- 
agates outward from i?phot(0- Detailed balance implies 
that the number of hydrogen atoms that are reionized 
x{z)nH{z) can be determined from (see §2.1 in Oster- 
brock 1988) 



(1 - x{z))nii 



AnjjCrEdE — x{z)'^n^a{T) , (7) 



13.6 gV 
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where j-^ is the photon flux in units of cm ^ s ^ sr ^ 
erg-i, a{T) = 4.18 x IQ-^^iT/W^ cm^ s'^ is 

the recombination coefRcient (MapeUi & Ferrara (2005)) 
and (Te ~ 6 X 10^^^ cm~^ is the photo-ionization cross- 
section of hydrogen atoms near threshold. 

The flux can be related to A'ion at a particular dis- 
tance from the QN. Since the thickness of the ionization 
boundary is much smaller than the Stromgren sphere, 
we can take the radius of the Stromgren sphere i?^ as 
the typical distance at which Eq. ([7]) applies. Using a 
typical parameter set Rs « lOOpc, i?phot(0=3x lO^^cm 
and Tpk=3 x lO^'^K, we then obtain the solution for x{z) 
from Eq. To determine the volume fraction f{z) 

of ionized hydrogen in the universe as a function of z, 
we use f{z) « i'{z)x{z) where the flUing factor of non- 
overlapping Stromgren spheres from QNe, I'iz), is deter- 
mined from (Barkana & Loeb (2001) 

du{z) ^ G.LTQNiV//B UF,oi{z) \ _ a(r)CngV(z) 
dz ~ 0.76 \ dz ) a^{l + z)H{z) 

In this equation, Nj denotes the number of ionizing 

photons released in a QN per baryon in the star, is 
the present-day density of neutral hydrogen and a(i) is 
the scale factor defined such that a(G) = 1. The factor 
of 0.76 in the denominator is the primordial mass frac- 
tion of H (we adopt 0.24 for He; see ^43]) . Fco\{z) is the 
collapse fraction at high redshift which is estimated sim- 
ply as the mass fraction in halos above the atomic cooling 
threshold in the Press-Schechter model (Press&Schechter 
1974). An order of magnitude estimate, assuming 10^^ 
baryons in a 1.4Mq neutron star gives iVj/^ = 1000 and 
yields 0.1a;QN-^^//_B = 10, which we keep fixed for all the 
curves in the figures below. 

The most important factor in determining the ioniza- 
tion state is clumping factor C. We choose C= 1,10, 30 in 
displaying results in Fig. [TJ These values are typical of 
the reionization era and are supported by more detailed 
modelling of the evolution of the inhomogeneities of the 
IGM (Furlanetto & Oh (2008a&b), Trac & Cen (2007), 
Miralda-Escude et al. (2000)). 

Our approximations also assume that a particular ion- 
ization front can be associated to a unique source, which 
is valid until the late stages of reionization when these 
fronts overlap. Since we do not include such effects, we 
simply set v(z)—\ when v{z) crosses 1, leading to the 
artefact of discontinuities at low redshift in Fig. [TJ 

The vertical dashed line shows the observational lower 
limit on the reionization redshift (Fan et al. 2003) z=5.8. 
For C < 10 (moderate recombination rate), it is clear 
that QNe can make a large contribution to reionization. 

We can now plot the results for the optical depth 
as a function of reionizing source turn-on sX, z = zq - 
shown in Fig. [2l These plots are determined from Eq. ([6]) 
for C = 1,10,30 with /re=10. The best-fit value for 
from ACDM cosmology (Dunkley et al. (2009)), assum- 
ing prompt reionization a,t z = zq, along with la limits 
are also shown. 

The contribution of QNe to the optical depth depends 
sensitively on the clumping factor®, but can be significant 
for C ^ 10, while for C ~ 30 or higher, the contribution 

* We have varied the clumping factor from 5-15 and find that cor- 




Z 

Fig. 1. — The fraction of ionized volume in the universe f{z) from 
Quark-Nova events for various clumping factors C=l,10,30. The 
vertical dashed line denotes the observed lower limit of reionization 
at a redshift of 2 = 5.8. 

is small. If C 1, including QNe leads to too large an 
optical depth in comparison to the 5-year WMAP data. 
Other reionizing sources such as Pop III stars have been 
studied for consistency with the WMAP optical depth, 
with most reasonable models of their evolution suggest- 
ing a "shortfall" in the optical depth (tc ^ 0.05) due 
to feedback effects that lead to self-regulation (Sokasian 
et al 2004) , even in the case of high escape efficiency of 
the ionizing photons. Based on the results displayed in 
Figsd] and [21 we suggest the intriguing possibility that 
QNe can help make up this shortfall. 
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Fig. 2. — The Thompson optical depth Te\z) from eq. JSJl for 
various clumping factors C=l,10,30. 



4.3. Helium reionization 

Helium reionization occured more recently, as sug- 
gested by a rapidly varying Hell optical depth at z « 2.9 
with a small line of sight variation Az k, 0.1 (Reimers 
et al. 2005,2006). Observations of the Hell Lya forest 
(Fardal et al. 1998) along lines of sight to bright quasars 

respondingly, varies from « 0.13 to 0.066, assuming prompt 
reionization at z ~ 20. Thus, we do find sensitivity to the clumping 
factor, but this senstivity diminishes at large values of the clump- 
ing factor. A similar observation is made in relation to reionization 
simulations in inhomegeneous IGM (Miralda-Escude, et al. 2000), 
where the effect of increasing overdensity saturates since late-time 
reionization (which is the era where the optical depth builds up) 
proceeds preferentially in underdense regions. 
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at z « 3 as well as empirical modelling of quasar lu- 
minosities and distributions (Furlanetto & Oh 2008a&b) 
hint that quasars are good candiates for Hell reioniza- 
tion. However, the latter work also cautions against a 
simple picture in which quasars are the only Hell ionizing 
source. This is because even a small fraction of HeH can 
absorb the redshifted Lya photons, so that a modest con- 
tribution from a widespread low-intensity source above 
the HeH edge can also explain the apparent patchy na- 
ture of Hell reionization at z > 3. Naturally, this raises 
the question: can QNe play a role in HeliumH reioniza- 
tion? 

First, we consider the case for Hel, which has an ioniza- 
tion potential of 24.6eV. We find from Eq.® that only 
5% of the photons that ionize HI can also ionize Hel; 
however there are far fewer Hel atoms nuc/nn ~ 0.08. 
Since the recombination rates of HI and Hel differ only 
slightly, the ionization front for Hel may end up closely 
following that of HI. If we assume that Helium reion- 
ization does not affect HI reionization, then for Tpk ~ 
3 X lO^K, we find that the radius ratio of Stromgren 
spheres R^^^yRf^ w 1. 

In reality, the two ionization fronts are coupled since 
photons above 26.4eV can ionize both H and He. While 
we have not performed a detailed solution of the cou- 
pled fronts here, the results of such an exercise can be 
mocked up by a parameter y = nHiaj/(HI)/(nHiai/(HI) + 
nHciaiy(Hel)) which is the fraction of photons with en- 
ergy E — hv > 24.6eV that are used up to ionize 
Hydrogen {n is the number density and is the pho- 
toionization cross-section). Although is strongly en- 
ergy dependent, the ratio a!y(HeI)/a^(HI) is about 8 for 
hv > 24.6eV. With nHoi/nm ~ 0.08, we find y ~ 0.65. 
This implies that {1 — y) ~ 0.35 is the fraction of such 
photons that ionize Helium. The ratio of HII to Hell 
Stromgren radii depends on y as well as on the emis- 
sivity of the source (and recombination coefficients) and 
the ratio R^^"/R^^^ « 0.35 for Tpk « 3 x lO^K. For a 
higher Tpk w 4 x 10'*K, this ratio is almost 1. It follows 
that the ratio of volume fractions V"^" /Vg" « 0.04 
for Tpk « 3 X lO-^K and « 1 for Tpk « 4 x lO^K. Includ- 
ing other smaller effects in Eq.(8) of our paper, such as 
the difference in recombination rates between H and He, 
and the number of ionizing photons/baryon for H and 
He, we find that the volume fraction of ionized Helium 
can vary dramatically from about 2% to 100% that of 
ionized Hydrogen. With such a strong dependence on 
the temperature of the ionizing source, it is not possible 
to make a robust claim on the efficacy of the Quark-nova 
on Hel reionization. 

For Hell ionization, with a large ionization threshold 
of 54.6 eV, the emissivity of QNe implies that only about 
1 in 10^ photons can reionize Hell. Employing the front- 
decoupling approximation as for Hel ionization, we find 
that i?f'="V-Rs" ^ 0-015 and that even for the largest 
?pk=4 X lO^K, approximately 5-10% of Hell is reionized 
by QNe by z ~ 3, depending on C ^. Thus, we can con- 
servatively say that QNe cannot be major sources of He- 
lium reionization but they may constitute a low-intensity 
source compared to bright Quasars. 

^ This number will be lower once coupling between the HII, Hell 
and HcIII fronts is introduced 



5. QUARK-NOVAE AND IGM ENRICHMENT 

The large value of Tg implies that the first stars ap- 
peared as early as z ~ 20, affecting the chemical evo- 
lution of the Galaxy and subsequent star-formation. A 
fingerprint of early chemical abundances has been pre- 
served in the extremely metal-poor stars ([Fe/H] is used 
here as a proxy for the metallicity) such as HE-0107- 
5240 with [Fc/H]=-5.3 (Christlieb ct al. 2004, Bessell et 
al. 2004) and CS 22949-037, with [Fe/H]=-4 (Depagne 
et al. 2002, Israelian et al. 2004) which constrains corre- 
lations between nucleosynthetic yields and the reionizing 
photon flux, since both are affected (the former some- 
what weakly) by the star's metallicity. Previous works 
(Daigne et al. 2004; Tumfinson al. 2004) concluded that 
a top-heavy (40-100Mq) IMF with a lifetime of 50-100 
Myr is required to simultaenously satisfy constraints set 
by the nucleosynthetic pattern in EMPs and early reion- 
ization. This requirement can change once the contribu- 
tion of QNe to reionization is included. At present, we 
lack a comprehensive model of chemical evolution that 
takes into account the contribution from QNe. However, 
we can assess their importance in a global sense. 

If most of the metals produced in supernovae explo- 
sions of the first stars are expelled into the IGM, Ricotti 
& Ostriker (2004) have estimated that the IGM metal- 
licity is given by 

ZiGM ^{l- /BH)3g X 10-3^0 (^) (9) 

where /bh is the collapse-fraction of massive stars into 
black holes, 0.3 < g < 2 depends weakly on the metal- 
licity, and Zq is the solar metallicity. This assumes that 
the first stars are the sole reionizing sources, in which 
case, one needs to invoke a large black hole collapse frac- 
tion ^ 0.3 in order to satisfy the observational constraint 
on the metallicity of the early IGM Zjqm — 2 x 10~'^Zq 
(Schaye et al.(2003)). Such a large collapse fraction re- 
quires a top-heavy IMF. 

QNe can change this scenario dramatically. As already 
shown, QNe can generate enough photons to provide a 
large optical depth, particularly if the clumping factor is 
small. Consistency with the optical depth then implies 
that the first stars contribute much less to the UV photon 
flux, so that Te from the first stars can be much less than 
0.1. From Eq.®, it is apparent that even with a small 
collapse fraction into black holes and a normal IMF, one 
may be able to satisfy the observational constraint on 
the IGM metallicity. It may not be necessary to invoke 
a large population of short-lived massive stars very early 
on and the IMF of Pop HI stars could actually be more 
like the present-day IMF. We also note that since QNe 
only produce elements above A ~ 130, they essentially 
do not contribute to metallicity, which is associated with 
the production of much lighter elements such as C and 
O. Early IGM metallicity and reionization are decoupled 
in the QN scenario. 

However, there is a link between early r-process abun- 
dance and reionization. The amount of r-process mate- 
rial ejected in a QN, as estimated from the energetics of 
the underlying phase transition, is about 1O~'*M0. The 
frequency of such events during the reionization era is 
/dsQN ~ -01 yr~^ Galaxy"^, so that the total r-process 
material ejected by z ~ 6 (10^ yrs) is about 10'^ M©. 
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Since the total r-process abundance in our Galaxy at 
present is constrained to be ^ IO^A'/q (Wallerstcin et 
al. (1997)), we conclude that about 10% of the total r- 
process elements at ^ > 130 was produced early on in 
QNe, with the rest produced in Type-II supernovae or 
binary neutron star mergers, as QNe became much less 
frequent. In this way, QNe also provide a prompt and 
local initial enrichment of heavy elements A > 130 that 
is seen in some of the oldest stars found to date. ^° We 
plan to study the implications for the observed scatter of 
r-process abundance in metal-poor stars in a subsequent 
study with a chemical evolution model. 

6. CONCLUSIONS AND DISCUSSION 

We have examined the role of Quark-Novae in reion- 
izing the universe and contributing to its metal enrich- 
ment. In particular, we have suggested that: 

(i) the most massive stars (M > AOMq) in a conven- 
tional or slightly top-heavy IMF collapse to black-holes 
with a possible (small) contribution to reionization from 
accretion but not to metallicity; 

(ii) the reionization is driven by intermediate- mass Pop 
III stars, whose higher mass members {25Mq < M < 
4OM0) end up as QNe, providing the bulk of reionzing 
photons and enriching their environment in elements be- 
yond A ~ 130; and 

(iii) low mass Pop III stars {8Mq < M < 25 Mq) end 
up as type II SNe. A long-lived and metal-poor popula- 
tion of low-mass stars begins to emerge at the end of the 
reionization epoch (Greif et al. (2008)). In our scenario, 
the dying out of the first heavy stars coincides with a 
peak in the QN rate and therefore a peak in ionizing 
radiation. 

Our calculations have shown that the photon flux pro- 
duced in dual-shock Quark-Novae can be large enough to 
overcome recombination in the reionization epoch. Com- 
plete reionization by 2: ~ 6 is achieved if the clumping 
factor of matter is small C ~ 0(1) but the optical depth 
is then ^ 0.17, almost twice that of the latest WMAP 
measurement . We have made estimates for the evolution 
of the ionized fraction of the IGM in a simple model of 
non-overlapping outward-propagating ionization spheres. 
Our results imply that the optical depth from Thomson 
rescattering of CMB photons can be close to the value 
measured by WMAP, if the clumping factor C ^ 10. In 
this case, QNe provide about 60% of the reionizing pho- 
tons. Alternatively, if C ~ 0(30) or larger, QNe play 
a minor role in reionization compared to the first stars. 
QNe can be important for HI and Hel reionization, but 
their spectra do not contain sufficient high-energy pho- 
tons to reionize Hell in any significant amount. Of cen- 
tral observational interest is the production of r-process 
elements beyond A ~ 130 in QNe, a feature which dis- 
tinguishes QNe from SNe (the latter require larger en- 
tropy than predicted in current simulations to produce 
r-process elements beyond A ~ 130). 

There are several novel aspects of the dsQN which must 
be addressed in future. The most important is the final 
composition of the Supernova ejecta from the first explo- 
sion which is subsequently shocked by r-process rich QN 

The large ratios of [C/Fc],[0/Fc],[N/Fc] seen in the EMPs 
is believed to be related to lower mass stars, 20 < M/Mq < 40 
and would not be affected by QN events, which require heavier 
progenitors and only produce elements beyond A ~ 130. 



ejecta. This is a critical input to chemical evolution stud- 
ies which can assess the effect of QNe on early r-process 
abundance patterns observed in metal-poor stars. In ad- 
dition, the effect of QNe on the extra-galactic radiation 
field as well as its potentially large impact on the dis- 
sociation of molecular Hydrogen, which is a "negative 
feedback" to the formation of massive stars, remains to 
be analyzed. Further details relevant to the reionization 
era, such as the metallicity, the evolution of the clump- 
ing factor, the star formation efficiency in Galaxies, and 
the early supernovae rate will have some quantitative im- 
pact on our results. These quantities are poorly known at 
present and present a source of uncertainty in any model 
of reionization. 

We close by predicting some important distinguishing 
features of QNe that are amenable to cosmological ob- 
servations: 

i) Quark-Novae and nucleosynthesis: QNe are a novel 
nucleosynthetic site for the r-process. We expect the QN 
ejecta to achieve 7-ray transparency sooner than Super- 
nova ejecta since QN progenitors (i.e., neutron stars) lack 
extended atmospheres. Thus r-process-only nuclei with 
7-decay lifetimes of the order of years (such as ^^^Cs, 
i44ce, i55Eu and ^^'^Os) can be used as tags for the QN 
(Jaikumar et al.(2007)), differentiating them from pair- 
instability SNe (due to the lack of any neutron excess) 
or core-collapse SNe (lower neutron excess than the QN 
ejecta). This could be observed by Gamma-ray satellites 
such as INTEGRAL. 

ii) Quark- Novae and high redshift gamma-ray bursts: 
It has been argued that QNe provides an additional 
energy reservoir in the context of Gamma Ray Bursts 
(GRBs). The QN leads to a three stage model for long 
GRBs, involving a neutron star phase, followed by a 
Quark- Star (QS) phase and a plausible third stage that 
occurs when the QS accretes enough material to become 
a black hole. As shown in Ouyed et al. (2007) and Staff et 
al. (2007) by including the QS phase, one can account for 
both energy and extended duration of the prompt emis- 
sion, X-ray flaring and the flattening observed in GRBs 
light curves. Our flndings in this paper link QNe to the 
reionization epoch. The connection between QNe and 
GRBs suggested above would imply that GRBs should 
be observed as far back as the epoch of reionization. 

To further pursue this connection between early QN 
and high redshift gamma ray sources, we note that the 
recent detection of GRB 080913 with Swift (Gehrels et 
al. 2004) at redshift 6.7 makes it the highest redshift 
GRB to date - more distant than the highest-redshift 
QSO (Fynbo et al. 2008). At ^ = 6.7 the burst occurred 
when the universe was less than a Gyr old when a high 
fraction of massive stars is expected to be of Pop III. It 
is thus possible that this GRB is a member of the long- 
soft GRBs produced by the collapse of the massive Pop 
III star (Greiner et al. 2009). However, making GRBs 
by the collapsar mechanism from Pop III has been ques- 
tioned since simulations suggest that the exploding stars 
will retain its hydrogen and helium rich outer layers (e.g. 
Lawlor et al. 2008 and references therein). This violates 
the main requirement in the collapsar scenario that the 
star loses its extended envelope to enable the relativis- 
tic jet to punch through the compact core on timescale 
for long-soft GRBs (Fryer et al. 2001). Pop III binary 
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evolution can both help remove the hydrogen envelope 
and spin up binary components. However, recent sim- 
ulations question the efficiency in producing GRBs in 
this process (e.g. Belczynski et al. 2007 and references 
therein). Furthermore, classifying it as a long duration 
GRB begs the question, in the coUapsar scenario, of how 
can a massive star produce a burst as short as 1 sec- 
ond? If instead GRB 080913 belongs to the short-hard 
class (Pal'shin et al. 2008; Xu 2008), then a neutron 
star-black hole merger is favored over a double neutron 
star merger, with the Blandford-Znajek process at play 
(Perez- Ramirez et al. 2008). 

The present detection rate of GRBs at 2; > 5 is about 
what is predicted on the basis of the star formation rate 
(Jakobsson et al. 2005). In addition, with hints of first 
stars having formed as early as z > 20 (Kogut et al. 
2003), GRBs are believed to exist as early as z ^ 15-20. 
The lack of many high- 2; GRBs can be explained in our 
model, as a consequence of clustering of GRB events near 
the end of the cosmic reionization era. We recall that in 
our model, the dying out of the first stars (at z ~ 7-8 
when adopting a normal IMF) coincides with a peak in 
the QN rate. This peak in the QNe rate (i.e. a peak in 
the GRB rate at 2; ~ 7-8 in our model) offers a possible 



explanation for the sparsity of GRBs out to the highest 
redshifts (z > 10) despite the immense luminosity of both 
the prompt gamma-ray emission and X-ray and optical 
afterglows of GRBs, and the current technology (i.e. the 
high detection rate delivered by Swift). While a merger 
origin of GRB080913 cannot be definitively ruled out, 
we suggest that this burst could instead be the first ever 
detection of a QN near the end of the cosmic reionization 
era. 
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APPENDIX 
A; POLLUTION BY CFL STRANGELETS 

By assuming that strange quark matter is absolutely stable, the universe might be polluted by strangelets (e.g. 
ejected during coalescence of quark stars in binary systems; e.g. Madsen 2005). With a NS-NS binary collision 
rate of 10"** — 10~^yr~^Galaxy~^ (Belczynski et al., 2002), we expect the rate of QS-NS collisions to be at most 
10~^ — 10~''yr~-^GaIaxy^-^. This will be somewhat reduced by the fact that, just as about 50% of binary systems are 
disrupted by a supernova, binary neutron stars can become gravitationally unbound in a quark-nova. The amount 
of ejected mass will also be smaller since quark matter is stiffer than neutron matter. Following Madsen (2005), we 
have estimated the flux of CFL strangelets in the Earth's vicinity, accounting for factors such as the geomagnetic 
cutoff on the rigidity of the strangelet, and the confinement time of typical CFL strangelets in our Galaxy. We 
find this flux (assuming strangelet number distribution to peak in the baryon number range A ^ 10^ — 10"^) to 
be ^ 10^ — 10'*m~^yr~^sr~^. This flux is too low to be tested conclusively by terrestrial experiments (although 
for particular values of the strangelet charge, 1 or 2 candidate events have been identified; see Finch 2006). The 
satellite-based AMS-02 detector launched recently has a threshold senstivity that is sufficient to detect the expected 
flux of CFL strangelets from quark star collisions across a wide range in A; these results are expected within a few 
years. Lunar searches, which benefit from lack of geological mixing and no magnetic field deflection, are also just 
approaching the required flux senstivity. Thus, at present, the strange quark matter hypothesis is not inconsistent 
with observations, even taking into account the binary QS-NS collisions and consequent pollution by CFL strangelets. 
The most recent work about strange star binary mergers by Bauswein et. al. (2008) found that for high values of the 
MIT bag constant, strange stars could co-exist with ordinary neutron stars as they are not converted by the capture of 
cosmic ray strangelets. Combining their simulations with recent estimates of stellar binary populations, Bauswein et 
al. (2008) conclude that an unambiguous detection of an ordinary neutron star would not rule out the strange matter 
hypothesis. 

B: QUARK-NOVAE PROGENITORS 

The fit to the observed light curve of SN2006gy, Leahy & Ouyed (2008) assumes QNe progenitor mass in the (40- 
60)Mq range. However, one can employ the parameter degeneracy in that fit to examine the dual-shock scenario with 
the more conservative mass range of (25-40)Mq which is more in line with the literature (e.g. Heger et al. ApJ 591 
2003; Nakazato et al. 2008) which suggets prompt BH formation above 40AfQ. It should be noted however that the 
effect of the fireball in the CFL phase (Ouyed et al. 2005) has not been taken into account in any of these simulations, 
which means the range 25-4OM0 could be somewhat underestimated. 

C: SHOOK EFFICIENCY IN DUAL-SHOCK QUARK-NOVAE 

According to Ouyed et al. (2007), the shock efficiency varies as p^j^^ with the mean SN envelope density given by 
Penv oc Menv./-Renv choose 4OM0 instead of 6OM0 for the progenitor of the SN, and demand the same efficiency, 

we find that the collision radius should be R40 = (40/60)^/^ x Rqq ^ 0.876 x i?6o and the delay time tdeiay,40 = 
0.876ideiay,60 ~ 0.876 X 15 days ~ 13 days. Similarly, Tpk,40 = (i?4o/i?6o)^/^rpk,60 = (0.876) V2Tpk,6o ~ 0.94Tpk,60- 



Quark-Novae, Reionization, and Early R-Process Elements 



9 



With these changes, we find a reduction factor of 0.6 in the total number of ionizing photons. This does not change our 
order of magnitude arguments on the ionization efficiency and we had in any case adopted a very conservative estimate 
for the number of ionizing photons initially by neglecting photons emitted after the ejecta becomes transparent. 
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